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NOTE ON THE USE OF GROUP THEORY IN 
ELEMENTARY TRIGONOMETRY 


By G. A. MILLER 


1. THE applications of group theory in the development of some of the 
fundamental parts of trigonometry have received considerable attention :* while 
some of these admit of interesting amplifications yet most of them relate to 
spherical trigonometry and are too complex to be immediately useful to the 


average teacher of elementary trigonometry.. The present note relates to 


developments which are generally treated near the beginning of a course in 
trigonometry. While the mode of treatment which is outlined may not be 
suitable for the beginner, it will doubtless enable the teacher to see many 
points in a new light and to observe relations which had escaped notice. 

The only group properties which are employed are those which belong 


to the cyclic group and the dihedral rotation group. 


these groups will make the reading easier, yet it is not essential, in view of 


the elementary character of the considerations. In fact, these trigonometric 


developments furnish one of the easiest means of acquiring a clear notion 
of the groups in question. This paper has very close contact with the 
writer's paper entitled “A new chapter in trigonometry”t and includes the 


results developed in the first part of that paper. 


* KE. Study, Abhandlungen der Sachsischen Gesellschaft der Wissenschaften, vol. 20 (1893), p. 
83; F. Meyer, Crelle, vol. 115 (1895), p. 209; G. Chisholm, Dissertation. Gottingen, 1895; 
O. Pand, Hamburger math. Gesellschaft, vol. 3 (1897), p. 290; G. A. Miller, Quarterly Journal 
of Mathematics, vol. 37 (1906), p. 226. 

tQuarterly Journal of Mathematics, loc. cit. 
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We assume as known that 


sin(90° — x) = cos 2, cos(90° — x) =sin 2, 





sin(45° — x) = 4V2(cosx—sin x), cos(45°— 2x) = §$ V2 (cos «+ sin x). 


It is proposed to derive from these formulas, with the help of some 
elementary group theory, expressions for the functions of fourteen other angles 
in terms of those of 2. 








The operation of deriving the angle 45° — x from x is of order 2, since 
45° — (45° — x) =x. Representing this operation by ¢ and that of taking 
the complement of any angle by c, we have @=c? = 1. The group }f, c{ 
generated by / and ¢ is therefore of the dihedral rotation type and its order is 
twice the order of fc.* As ¢c signifies that an angle is subtracted from 45° and 
the complement of the resulting angle is obtained, it follows that fe » = x + 45°, 
and therefore fc is of order 8. Hence }¢, ¢/ is of order 16. 

If tec = ¢,, the operations of this group may be written as follows: 
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All the operations of the second row are of order 2 while those of the first 
constitute the cyclic group of order 8. The operations of the first four 
columns constitute the group of order 8 known as the octic group or the 
group of the square. This subgroup is generated by c and Gc = s, where « 
represents the operation of taking the supplement of an angle. If we apply 
these 16 operations to the resulting angles in order (since angles are taken 
modulo 27r) are as follows : 
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x e+ HP x + 1RO 2+ 270° 2 4+ 45° 2+ 185° 2+ 225° x2 + 315° 


90°—-x —2x 270° —x 180°—r 45°—2 315°—2 225°-2 135°-—2 


If these 16 operations were applied to any other one of these angles the same 


16 angles would result (in different arrangement), since the sixteen opera- 
tions constitute a group. 





* Bulletin of the American Mathematical Society, vol. 7 (1901), p. 424. 
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2. The trigonometric functions of sixteen angles. What pre- 
cedes has been preliminary to expressing the trigonometric functions of each of 
these sixteen angles in terms of those of « assuming only that the equations 
giving the functions of 45° — x and of 90° — x in terms of those of x are 
known. It will be sufficient to do this for the sine and the cosine since the 
other functions are readily derived from these two. We begin with # + 90° 
and for clearness give all the details in this instance. 


sin(x + 90°) = sin Gx = sinlctcx = cos etx = 4 V2 (cos lex + sin tex) 
= 3 ¥2(sin tx + cos tr) = 3 (cosx— sin x + cos x + sin x)= cos 2, 
cos(x + 90°)= cos fix = sin tctx =$Y2(cos tex = sin tex) = 3) 2(sin da — cos tx) 
= 4(cos x — sin x — cos x — sin x) = — sin z. 


The functions of the angles corresponding to ¢; and ¢; may be directly 
derived from the preceding results as follows : 


sin (x + 180°) = sin # = sin Gfx = cos fj = — sin 2, 
cos(x + 180°) = cos x = cos Gf«# = — sin Gx = — cos 2; 
sin (x + 270°) = sin fj = sin Gtjx = — sin fix = — cos 2, 
cos(x + 270°) = cos fi.“ = cos (Gx = — sin x= sine. 


The functions of the four angles corresponding to c, Gc, fe, Ge can be 
directly written down by means of our knowledge of the effect of the funda- 
mental operation c and the results which have just been obtained. The work 
may be done as follows : 


sin(90° — x)= sin cx = cos 2, 

cos(90° — x)= coscx = sin x; 
sin(— x) = sin &\cx = cos fx = — sin 2, 
cos(— x)= cos ficx = sin §x= cos az; 
sin (270° — x) = sin fjcx = cos tix = — cos x, 
cos(270° — x) = cos Gjcx = sin (7 = — sin x; 
sin (180° — x) = sin fer =cosGx = sina, 
cos(180° — x) = cos fcx = sin fx = — cos x. 
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The functions of the three angles corresponding to 4, 4, 4 may be ex- 


pressed as follows : 


sin (2 + 135°)=sin G@r=sinhhxe= coshe= sin fx 4V¥2(cosa—sin2z), 
cos (x + 135°)= cos fix = cos Off = — singe =— cos tr = — 4 ¥2(cosx + sin x) ; 
sin (vv + 225°)=sin Ga =sin Ofje = — sin hw = — cos tr = — $ ¥2(cosx + sin), 
cos (.¢ + 225°)= cos Gv = cos Afr coshw=—sinte= 492(sine—cose) ; 


sin (vw + 315°)=sin Gv =sin (fe = —costx=—sintv= 4$/2(sinx—cos.), 


‘ 


cos(#+315°)=cos fx=costfx= singw= costr= 3 V2(cosx+sinz). 


As the functions of the angles which correspond to the remaining opera- 
tions of our group may be obtained by first passing to the co-function and then 
proceeding as above it seems unnecessary to give them here. The thing which 
we desire to emphasize is the method of deriving these functions and that a 
knowledge of this group of order 16 gives us a clear insight into the number 
of possible similar ways of attaining the same results. Since this group is 
generated by 47 other pairs of operations we could have chosen the two fun- 
damental angles in 48 distinct ways * and the derivation of the trigonometric 
functions of the remaining fourteen angles could have proceeded along the same 
lines. This clear comprehension of the various roads which lead to the same 
goal combined with the freedom to choose from such a variety of almost 
equally good roads constitutes one of the greatest sources of pleasure to the 
appreciative traveler in a mathematical field. 

If it had been desired to find only the functions of the angles corre- 
sponding to the octic group we could have used the angles which correspond 
to two of its generators as fundamental instead of those corresponding to c 
and ¢. This was done in the article published in the Quarterly Journal to 
which we have already referred. The method pursued here presents the mat- 
ter under a more comprehensive aspect and exhibits more extensive relations. 
There is however no limit to the generalizations along this line. For instance, 
if the operation ¢ were replaced by an operation 7, where rz indicates the result 
of subtracting the angle z from 30°; rz = 30° — «, we should have to deal with 
a dihedral rotation group of order 24 which would again contain the octic group 


*The number of ways of choosing two generators of the dihedral rotation group of 
order gy is 49 ¢ (4g), o(n) being the totient, or indicator, of n. 
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asa subgroup. The considerations with respect to this group are practically 
the same ag those given above. 


3. The trigonometric functions of forty-eight angles. Since 
rtx = x + 15°, its order is 24 and jr, ¢! is the dihedral rotation group of 
order 48. The angles which correspond to this group are x increased by the 
twenty-four multiples of 15° together with their negatives. The functions of 
these 48 angles can therefore be directly expressed in terms of those of x pro- 
vided we know those of the angles corresponding to generators of jr, t}. As 
we have assumed a knowledge of those which correspond to ¢ it is only neces- 
sary to add those which correspond to r; i. e. ‘ 
sin (30° —x) = 3(cosx—-¥3 sin x), cos (30° — x) = $(¥3 cosx+sin 2). 

Hence, 

sin(x + 15°) = sin rtx = 4 ¥2 (cos rx — sin rz) 
= } ¥2(V3 cos x + sinz — cos x + V3 sin x) 
= ¢ V2} (V3 + 1) sin + (V3 —1)cosz{, 
cos(#+15°) = cos rtx = 4 V2 (cos rz + sin rz) 
=4 V¥2(¥3 cos x+sinxz+cosx — V3 sin x) 
= $25 (V3 + 1)cos x —(¥3 — 1) sin xt. 

The functions of the remaining 47 angles may be readily found from these 
results if we observe that }r, t} = jrt, c{, and hence rt = ¢, andc may be 
used as the fundamental operations. To obtain the functions of 75° — x= rtcx 
it is only necessary to observe that sin rfcx = cos rtx and cos rtcx = sin riz, 
from which it follows that 
sin (75° — x) = sin ¢,cx = cos (2 = $ V2} (V3 + 1)cos a — (¥3 — 1)sin x}, 
cos(75° — x)= cos f,ex = sin t,x = $ V2! (V3 + 1)sin x + (V3 — 1)cos x}. 

Further, 
sin (x + 30°)= sin Gx = 4 V2} (V3+ 1) sin 42 +(V3 — 1)cos t,x} 
=$}(4+4+2 V3)sin x + 2cosx+2cosx—(4—2 V3)sin «(= 4(cos 2+ V3sin x), 
cos(x + 30°)= cos Gx = V2} (V3 + 1) cos &x — (V3 — 1) sin tz} 
=} (4+ 2V3)cosx— 2sin x— 2sinx —(4— 2V3)cosa{ = $(V3 cos x—sin 2) ; 
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sin (2 + 60°) =sintz = 4(cos Gx + V3sin Gx) , 
= 4(V3 cos x—sinz + V3cosx + 3sin x)= $ V3 cos x + sinz), 
cos(x + 60°) = cos 4x = $(V3 cos Gx — sin Gx) 


= }(3cos x— V3sin x—cosxz — V3 sin x)= 4(cosx—V3 sinz). 


The preceding examples appear ample to exhibit the method of work and 
the numerous relations which it sets forth. The two generating operations of 
‘r, ¢{ could be selected in 36-8 = 288 different ways so that there is an almost 
endless variety of ways from which to choose in deriving the functions of these 
48 angles along this line. The method aids in deducing a great deal from 
a few known facts but is of no value without such knowledge. It has a theo- 
retic and clarifying influence but does not seem to fit in a brief practical course’ 
in trigonometry. If the trigonometric functions of a series of angles is known, 
those of any multiple (positive or negative) of the highest common factor of 
all the possible differences between these angles may be found in this way.* 
In fact, we are dealing simply with the group of subtraction from several 
fixed numbers with respect to modulus 27, which is considered in the last 
article to which we have referred. The present article may serve as an illus- 
tration of the utility of such groups. 





* ANNALS OF MATHEMATICS, ser. 2, vol. 6 (1905), p. 41. 








THE ASYMPTOTIC LINES ON THE ANCHOR RING 
By Marion BALLANTYNE WHITE 


THE purpose of this paper is to determine the equations of the asymptotic 
lines on the surface of the anchor ring when its equations are expressed in para- 
meter form and, with the aid of elliptic functions, to examine their geometrical 
properties.* The first three sections are devoted to the equations of asymptotic 
lines in general, to their applications to the anchor ring, and to the ex- 
pression of the equations of the asymptotic lines in terms of the elliptic 
functions of Jacoby. After determining the general form of the asymptotic 
curves, the anchor rings for which they are closed will be discussed. 
The paper concludes with the computation of a numerical] example. 


1. Asymptotic Lines in General. Consider a surface whose equa- 
tions are 
xr=x(u,v), y=y(u,v), z=2(u, v), 


and a plane tangent to that surface at a point .V whose coordinates are (u,v). 
The distance of any point (« + 4, v + &) on the surface in the neighborhood 
of .V from this tangent plane is represented in Bianchi’s notationt by 


(1) §=3(Dh? + 2Dhk + D'F*) +0, 





*The only literature I have found on this subject is a brief summary of an article by 
T. Matoda, ‘‘Asymptotic lines on a circular ring,” Tokio Math. Soc., vol. 4, pp. 217-219, in 
the Jahrbuch ther die Fortschritte der Mathematik, vol. 21, p. 765, and also an article by Prof. 
Maschke, ‘‘Asymptotic lines on acircular ring,” Bull. Amer. Math. Soc., vol. 2, 1895, pp. 19-21, 
in which the equations of the lines are briefly reduced to a form involving Weierstrassian ellip- 
tic functions without any corresponding geometrical investigation. 

tL. Bianchi, Vorlesungen ther Differentialgeometrie, translation by M. Lukat (1896), 
pp. 61, 87, 104. 
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where 7 is an infinitesimal of the third order in 4 and &. The expressions for 


D, D’ and D” are: 
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where x,, x, are the partial derivatives of x with respect to u and v, and as 
usual 
(3) B=e+94+%, Parry ret YuYot Mure GF=HMyt Yo + % 
The surface in the neighborhood of © lies on one side of, or is cut by, 
the tangent plane according as 6 retains or changes its sign. This 6 will 
always be positive for sufficiently small values of / and & if the two factors of 
Di? + 2D'hk + De ave conjugate imaginary factors, i. e., if DD" — D? > 0; 
in this case the surface is said to have elliptic curvature at the point M. If 
DD" — D® < 0, & may take positive or negative values for small values of A 
and / and the curvature at W is said to be hyperbolic. Parabolic curvature 
occurs at points where DD” — D?® = 0. 
The curve in the uv-plane corresponding to the equation 


(4) 8 = 3(Dh? + 2D/hk + D2), 


where 6 is constant and 7 has been neglected, is approximately the image in 
the uv-plane of the section of the surface by a plane parallel to, and at a dis- 
tance 6 from the tangent plane. 

Two systems of curves on a surface are conjugate systems if the directions 
of their tangents at every point on the surface correspond to directions in the 
urv-plane which are conjugate with respect to the conic (4). From the theory 
of conic sections, the condition for this is* 





* See also: Bianchi, Vorlesungen ther Differentialgeometrie, p. 108. 
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(5) Ddu Su + D'(du bv + dv bu) + D'dv dv = 0, 


where d and 6 are used to distinguish between the increments along the two 
curves. In particular, if a curve on the surface is such that the direction of the 
tangent at any point is conjugate to itself, the curve is called an asymptotic 
line on the surface. It follows immediately from the last equation that the dif- 
ferential equation for the asymptotic lines is 


Dd? + 2D'du dv + D'dv? = 0. 


As this is a quadratic in dv and dv, the asymptotic lines evidently forma 
double system. If they are real, the equation must have real roots. It fol- 
lows therefore that the real asymptotic lines lie only on that part of the surface 
where DD” — D” < 0. 

To summarize: The asymptotic lines on the surface x = x(u, v), 
y=y(u,v), z=2(u, v) forma double system of curves whose differential 
equation is 
(6) Dadi? + 2D’ dudv + D' dv? = 0, 


D, D' and D" having the values given in equations (2). They are real 
only on that part of the surface where DD" — D® < 0. 


2. The Equations of the Asymptotic Lines on the Anchor 
Ring. Let us now apply the general formulas just obtained to the particular 
case of the anchor ring. The anchor ring to be considered is generated 
by the revolution of a circle with radius 4 and center at (a, 0, 0) about the 
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Fic. 1. FIG. 2. 


z-axis. Instead of using the equation in the form F(x, y, 7) = 0, we may ex- 
press 2, y and zas functions of two parameters rand v or of and v, whose 
meaning is indicated in the accompanying figures. 
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By referring to these figures, it is easy to see that the equations of the 
anchor ring in terms of u, v are: 
(7) x= (a—hecosu) cos v, 
y = (a — bcos #) sin v, 
z —hsinu, 
the parameter lines being the meridians and parallels on the surface. 
By means of the formulas (2) and (3), 2, FP, G@, D, D’ and D” may be 
calculated, with the following results :— 
E=?, F=0, G=(a—+heosu)?, VEG — F? =h(a — 6 cos u), 
D=-A, D=0, dD" = (a—h cos u) COs WU. 
It is evident that F’ must equal zero, since 7,, Y,,. 2, and 2, Y,, Z, are 
the direction cosines of the parameter lines, which intersect orthogonally. 


Substituting these values of 1), J, and D” in the equation (6) of the 
asymptotic lines, it becomes 


(8) — hbdi? + (a —h cos u) cos undy* = 0, 


from which we derive 








(9) dv = + Vb —_—"s 
Vcos “(a — beos 1) 

or 

(10) vot V/, [ —— a..2 ee 
Jy Veosu(a —heosu) 


The constant C’ may be determined so that the asymptotic line passes 
through an arbitrary point (%, 7) on the ring, provided (1%, 7%) lies in the 
portion of the surface where DD" — D? <0. Then 





(11) C= F v2 | a... ene 
Jy Veosu(a — heosu) 
and 
u ] “u,, 
Jo Voosu(a — heosn) » Ycosu(a — bcos u) 








1907] THE ASYMPTOTIC LINES ON THE ANCHOR RING 107 


By a rotation about the z-axis, the equation for any asymptotic line can 
therefore be finally reduced to one of the forms 


a du 
~— wore Vcos u(a — b cos u) 


In §1 we found that these lines can be real only when DD" — D? = 
— bcos u (a—b cos u) < 0. Since a > 4, we have a — 4 cos u > O and there- 








fore this inequality holds only when — 5 <u< a These are exactly the values 


for which the denominator in the integral of (13) is real. 
Our results are then as follows: 


1) The equations of any asymptotic lines on the anchor ring may, by a 
rotation about the z-axis, be reduced to one of the forms (13). 

2) The asymptotic lines are real only on the inner half of the ring, that 
~ » ;' ca 7 
is, in the region — 9 << > 

3. The Equations in terms of Elliptic Functions. The pur- 
pose of the present section is to express the equations of the asymptotic lines 
in terms of the elliptic functions. That the integral(13)is really an elliptic in- 
tegral, appears when wu is expressed in terms of r = a — 6 cos u (see figure 1). 
Since cos u = (a — r)/b, the integral (13) becomes 





v dr 
9 = ] eee . 
(38) ’ - = r)((a—r)?— 8} 


Hence r, and therefore also u, are elliptic functions of v. 
To reduce this integral to the normal form, Durége’s method* will be 


used. The integral to be reduced is in the form | 


r 





= where A = 1 and 


1 « J 


the four roots of 2? = 0 are 0, a — b, a, anda +. Weill now introduce a 
new variable &, so chosen that the values of & which correspond to 0, a — 6, 
a, and a + bare — 1/k, — 1, 1, and 1// respectively. If the four roots of the 


*H. Durége, Theorie der elliptischen Functionen (1887), p. 34. 
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equation in 7 are designated by 7, p, 9, and «, we have the following corre- 
spondence between the values of & and 7. 


= + 1,-—1, —1/k, 
- — 


where r =a+b, p=a, q=a-h, «=9. 


7, Ps Gs K, 


The following formulas are the ones given by Durége to effect the trans- 
formation to the normal form when the roots of #? = 0 are all real. The re- 
lation between r and € is given by® 


Pop gee 


() fae ten" 
where 
(17) pa../(P—7)(P — *), 


~=N(q—7)(q—*) 


The value of 4 is determined from 





1—k_ (p—7)(q—*) 
_— ltk ~ N(q—7)(p—«) 


The values of the old and new integrals are related by the equation 


webs _ dr oe 2Vh dé 
( *) Jar—pyr—or—-wr—e) Zz \VA( p—q)(7—«) V1 —-&)(1 — PR) 








or 
ly ? . 
(0) == ‘ a == ae = = ey 
VA(r—pyr—g)(r—m(r—h) VAlp—g)\r—K«) VI—Rsin*d 
where 
(21) — = sin ¢. 


On making the substitutions for the anchor ring from (15) in (16), (17), 
(18), and (20), we obtain the following results : 


25 bP(1— &) 
22 a. > ——— ——_$——___ -— 
ie "= “+ TE PO 6’ 
23 Nie a a 
— \ 2(a— hb)’ 








* Durége, loc. cit., p. 44, ff. 
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the minus sign being chosen, since it is desired that r and £ increase simul- 
taneously. Also 























(24) p= V24- lf 
V2a + Va —b 
and 
bk ff dg bk [*___ a 
25 '= +2 [_OR : poco Ss i 24/ ir 
(25) e=42V TT Javd—-Pd— Pe) a+b J :Vi— Paint 
Let 
y bk 1 
= 2/3 +b a 
Then in the notation of Durége,* equation (25) may be written 
1 f¢dd 
an nad (tt, 
a J_® Ad 
or 
dp °° do +d ‘ 
28 save [E+ az =/ ~z,+ 
sie w= if ag* ls ag~f ag * 


or, solving for ¢, 
(29) ¢?=am (t+av—F). 
The formula for calculating £ = sin ¢ ist 


snucnv dnv — snv cnu dnu 
1 — A*sn*usn?v 





sn(u—v)= 


’ 


u, v in this case having the values + av, A. After using the relationcn H= 0 
and sn A = 1, this gives 

_ _en(+av)jdn(+ av) enav 
(39) aie dn?(+ av) ~  dnav’ 








since cnav and dn av are even functions. 
It now remains to express u in terms of these elliptic functions. With 
the help of equation (22), since cos u = (a — r)/6 we obtain 


~P(l-£) 
1+&+P(1 - &)° 





(31) cos u = 





* Loc. cit., pp. 8, 17. 
¢t Durdge, loc. cit., pp. 106, 116. 
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Substituting (30) in (31) we find that the equations defining u in terms of v, 
which determine the asymptotic lines on the anchor ring, may be expressed in 
the form: 





{ P(dnav + enav) 
(4) cos = — aay — cnav — P(dnav + cnav)’ 
(32) 
(6) sinu= + \ (dnav— cn av)? —2 P(dn? av — cn* av) 
: dnav — cnav — P(dnav + cnav) 


where the values of P and a ure to he taken from equations (23) and (26). 


4. The Geometric Character of the Asymptotic Lines. It 
has already been shown that the asymptotic lines lie on the inner half of the 
anchor ring. Their form may be determined by examining equations (32) 
and (9) which give uv and dv dv respectively, and also the equation 


(33) dtu — sin u(24.cosu — a) 
dy? Qh 


obtained by taking the v-derivative of (9). 

In equation (52, 4), the + signs indicate that the asymptotic lines are 
related in pairs symmetrical with respect to the v-axis. Only the positive 
sign need be considered. The sign of ({) is not necessarily the same as that 
chosen in (32,7). The correspondence will be evident from the following 
discussion. 

The positive function uv detined by equations (32) is, in the first place, 
periodic in v since the elliptic functions have a period of 4A, where A is 
defined by 

"2d 


34 "= , 
(34) w= [xo 
The period in av is 4A, that is 


4h 
n(n 4 ) = EU). 
a 


Therefore to determine the graph of the line, we need only consider u as a 
function of v in the interval 0 Su = 4Aa. This interval can be still further 
restricted since 


en(2Hiw)=—enyw, dn(2Hiw) =dnw, 
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(= ) ‘3 
u +v) =u{— — v). 
a a 


The curve is thus symmetrical with respect to the ordinate v = 2/7/a, and it 
follows that its form is completely determined when the function u has been 
investigated in the interval 0 S v S 2Ff/a. 


and therefore 





Equations (32) and (9) show that w is a continuous function of v and has 
a continuous derivative. Moreover, du/dv can vanish only when cos u = 0, that 
is when cn av + dn av = 0, and it can be shown that this equation is true in 
the interval 0 = » S 2///a, only foray = 2H. The limits between which u 
varies can be found by putting » = 0 and v = 2A’/a successively in equations 
(32). When av = 0, since dn(0) = 1 and en(0) = 1, it follows that u = 0. 
For av = 2K, we have dn av = 1, cn av = — 1, and therefore u = 7/2. It is 
evident that for the positive function ~ defined by equation (32), the sign in 
equation (9) must be the positive one, and that uw constantly increases as v 
varies from 0 to 2//a. 

The equation (33) shows that the curve is always concave toward the 
v-axis in the interval 0 S v S 2/’/a, unless 24 >a. In any case there is an 
inflexion point at » = 0, and when 2 4 > a, there is one also at vu = cos—!(a/2d), 
at which the curvature changes from convex to concave toward the v-axis. 

From the foregoing discussion it is evident that the images in the uv 
plane of the asymptotic lines have the properties shown in figure 3. 











~ -”©6~—<CétY ys = se 
“iy, 
ae, o "Toole 


Fic. 8. 


On the anchor ring (see figure 4) the asymptotic lines oscillate back and 
forth across the inner equator between the two circles u = 1/2 and u = — 7/2. 
Any two asymptotic lines ave similar and can be made to coincide by a ro- 


tution about the axis of the ring. A section of one of the lines between 
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two successive contacts with the circle u=— 1/2 is symmetrical with re- 
spect to a plane through its point of contact with the circle u = 1/2 and the 
axis of the ring, and any such section may he made to coincide with any other 
by a rotation about the same axis. 


The period of the positive function u(v) defined by equations (32) was 
found to be 4A/a. But if the positive values of « in the interval 0 < v <4H/a 













are taken with the negative values in 44/a < v < 8A/a, and so on as indicated 
in figure 3, so that the line defined has a continuously turning tangent, the 
resulting function of » will evidently have the period 8/‘/a. 


5. Anchor Rings with Closed Asymptotic Lines. The problem 
here is to determine what form the anchor ring must take if the asymptotic 
lines are to be closed curves. Since they have a period of 8A/a, the condition 
to be satisfied is that some multiple of 8/7/a must be a divisor of a multiple of 
27, or stating it algebraically, 


8m 


a 


= 2rn, 
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where m and x are positive integers. From this equation, since A has the 
value (34), the condition becomes 


2 
an 7 
4m ah) V1 —Fsintd 





The right hand member depends on a and & and, as is evident from equations 
(24) and (26), may be considered a function of the ratio 4/a. Our first con- 
clusion then is : 

The value of the ratio h/a determines whether or not the asymptotic lines 
are closed ; i. e. the closure of the lines depends on the proportions and not on 
the size of the ring. 

The problem now is to determine whether there is any fixed ratio between 
a and } which will make the right-hand member of (35) equal to the constant 
mnj/im. We will first determine aand / in terms of where w = 4a. Using 
the values of 4 and a in (24) and (26). Then 


V2a—Va—h y2-Vl—w 





36 k= Pn a 

( ) V2a +vya-—b vi + Vi — ow 
- a+b lw +1 

(37) a=ty a = tyr’ 


If we designate by $(7), the furiction in the right-hand member of (35) ,the 
problem is to solve the equation, 


Tn 
ow) = >— . 


For this purpose we will study the graph of @(~), and since on the anchor 
ring @ > b, we may confine ourselves to the interval in which 0 < w<1. For 
w = 0, p(w) = 0, and when wv = 1, $(17) = ®, since’ = 1 and the integral 
in (35) is infinite. Therefore the curve passes through the origin and has 
w = 1 as an asymptote. 

The derivative ¢$'(7#) has the value 


, 1 da [{"? d 1/7? ksin*ddd dk 
(38) #(w)=— 3-5 fet ef ea de 
a 0 vyl—A*sin?d Fo (V1 —k*sin?g)” au 


2 du. 


i at i a a — 
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From (36) we derive 
dk v2 
~~ dw yl—w(y2+yl—w)® 


In (37) a is a function of both w and 4; therefore 


da Ca Ca ak 
(49) ae = Se * ER de 


l 


tf), fire roca 
4yh Fen +v) Now k¥l — w(¥2+ V1 —w)? 
Substituting these values from (3%) and (40) in equation (38), it appears 
that ¢'(7) has the value 








: [ 1, flu v2 a 
say k Ly3(1 + w) Vw kyl —w(V2+\ =f V1 — A*sin*¢ 


1 V2 r2 ksin*p dd 
ayi— w(y2 +4 eri (V1 — A sin?g)” 


— 





which is positive for 0 < ~ < 1, and therefore $(7) increases monotonically 
in that interval. 

The general form of the graph of @(1) is that of the curve (1) in 
figure 5. The straight line ¢ = ra 4m is represented by (2). Curves (1) 














FIG. 5, 


and (2) have but one point of intersection, which indicates that for every value 
of mand n that may be chosen, there is a fized ratio between a and b for which 
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the asymptotic lines will be closed: i.e., the proportions for a ring may be so 
chosen that the asymptotic lines close after any arbitrary even number 2m of 
intersections with the inner equator, and an arbitrary number n of windings 
about the z-axis. 


6. A Numerical Example. In this section it will be shown how 
the branch of the curve given in Figure 3, may be plotted from actual com- 
putation. The constants a = 18, 4 = 14 have been chosen so as to illustrate 
the case where there is a point of inflection for cos—! (a/2b). The values k = 0.5 
and a = 4/V14 are readily computed from equations (24) and (26). 

To find A, Landen’s first transformation* has been used. The object of 
this transformation is to increase the modulus & by successive steps, the upper 
limit of the integral in (34) being decreased incidentally until / is so near to 
unity that the integral has the approximate value 


The first transformation takes A’ into 


- 2 $y d¢, 
K = SS >=? 
l+h fy V1 — Ajsin*¢, 


where , and ¢, are determined by the equations 


—_— 2vk sin(2¢, — @)= & sin 
aha eo n( d, ¢) 5 d 


After two applications of the formula, A, is so near to unity that A’ may be 
written in the form 


Ih, 
k = V 7 log tan (2 + 1) : 


where - + + = 73° 41' 1.2”. Hence we tind A = 1.686. 








*W. E. Byerly, Jnteyral Calculus, p. 215. 
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The period of the asymptotic lines was found to be s#/a, which in the 
present problem therefore becomes 12.61. 
From the equations (32) and (9) of the asymptotic lines, values of u and 


corresponding to special values of v may be computed as shown in the 


following table : 








du 
? “ 
dr 
0 0 0.524 
Av ss ’ 
we 8.808 0.927 O.6415 
a 
2h 
— = a, Oe Tr 2 0 
a 











There is a point of inflection when uv = 0.872 radians. 
To find other points (vw, v7) on the curve, values of the integral 


= \ i 
? = | du, 
Jo Ycosu()} — 7 cosu) 


obtained from equation (13), may be computed by means of the trapezoidal 
rule for computing areas under a curve. 

This rule is: * “Add together the halves of the extreme ordinates and 
the whole of the intermediate ordinates, and multiply the result by the 
common interval.” To apply the rule, the interval 7/2 was divided into eigh- 
teen subdivisions of 0.087 each, the ordinate at each division point calculated 
and substituted in the formula just given. A table is appended which gives the 
resulting values of “and vy in radians; by means of which the curve in figure 
3, may be easily plotted for the interval 0 = v = 2A/a. 








* Williamson, Integral Calculus, p. 212. 
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For the projection of the asymptotic line on the xy plane, shown in 











FIG. 6. 


figure 6, the values for r were calculated from r = a — } cos u, and these are 


given in the third column of the table. 








TABLE 
u vu , 
000 .000 4.00 
087 .162 4.04 
174 .330 4.20 
.262 4381 4.46 
349 637 4.84 
.436 .788 5.30 
.524 934 5.80 
611 1.070 6.30 
698 1.206 7.23 
.785 1.350 8.10 
873 1.480 9.10 
960 1.620 9.80 
1.047 1.762 11.00 
1.134 1.804 12.10 
1.228 2.050 13.20 
1.309 2.203 14.30 
1.396 2.396 15.50 
1.483 2.631 16.70 
1.570 3.152 18.00 











UniversiTy oF Wisconsin, AuGustT, 1906. 
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ON CERTAIN THEOREMS OF MEAN VALUE FOR ANALYTIC 
FUNCTIONS OF A COMPLEX VARIABLE 


By D. R. Curtiss 


On account of the fundamental importance of the theorems of mean value 
for real functions of a real variable an investigation of the validity of the 
corresponding formulas for analytic functions of a complex variable may prove 
of interest, even if the results obtained are of no great value in themselves. 

Darboux has established theorems closely analogous to those expressed 
by the formulas : 


(1) F(z) — f(a) = (4 — a) f' [a+ 8 (xe —a)) 


=, - , [2 (|9@-43|<4) 
(2) [rw d(y)dy = f lu -8(x—a)] [o(u)dy 


but with the essential difference that in the second members of the above 
equations there appears a factor A whose absolute value is less than 1.* So far 
as I have been able to ascertain, no investigation of (1),(2), and allied formulas 
for real functions ofa real variable has been published which considers their 
validity without change of form when x, 7, @ and a take on complex values, 
and fand ¢ are analytic functions of their arguments. 

We may regard (1) and (2) as equations detining @ as an implicit fune- 
tion of x. When fand ¢ are real functions of a real variable subject to certain 
restrictions, at least one of the branches of this implicit function takes on areal 
value between 0 and 1 for each value of x in a given interval. The question 
arises whether equations (1) and (2) have still a solution @ of which at least 
one branch is always to be found in a certain limited region, when fand ¢ are 
functions of a complex variable analytic ina region / which contains the point 
a. We shall prove that there exists a neighborhood of « within which formulas 
(1) and (2), as well as others similar to these, remain true. In the case of 
real functions of a real variable something can be said as to the extent of this 


* Journal de Mathématiques, ser. 3, vol. 2 (1876), pp. 291, 294. Cf. Stolz, Differential- 
und Integralrechnung, part II, pp. 92-95. A recent publication by Brann — Beziehungen des 
Du Bois-Reymondschen Mittelwertsatzes zur Ovaltheorie—gives some interesting theorems 
which are not, however, entirely in terms of functions of a complex variable. 

(118) 
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region; in the present paper, where only analytic functions of a complex 
variable are considered, no positive results of this kind have been obtained. 


1, The Mean Value Theorem for Derivatives. The following 
theorem corresponds to formula (1) : 


THEOREM I. Let a be an interior point of a region B throughout which 
T(z) is analytic. Then there exists a neighborhood A of a such that for any 
value of z in A the equation 


F(6, z) = f(z) —y(a) — (2-4) f' [a+ @(z-—a)]=0 
is satisfied by a value of 0 which verifies the inequality 
6 —_ 4 < 4. 


To prove this theorem we expand /'(@, z) in a series proceeding accord- 
ing to powers of z—«. We thus obtain a development 


Mm: 


' F(0. = I hill L : 
(3) (6,z) = E - ai" (a)(z— a)", 
which converges for all values of z within a circle C whose center is a, pro- 
vided @ satisfies the inequality |@ — 4) < 4- 

If f(z) has the form az + p, the above theorem is evidently true since 
F’'(@, z) then vanishes identically. In all other cases, at least one of the suc- 


k 


cessive derivatives of f(z), beginning with the second, will not vanish at a. 
Accordingly we have 


Sa) =f"(a) = + =S"-9 (a) = 0, F(a) #0. 


We now make use of an auxiliary function ®(@, z), defined as follows : 


(6, z) = a z¢4, 
i @—«) 
1 Q@nr-1 ia 
(0,0) = (G- Goi): 


For this function we have 
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where "\ 1 is the real solution, lying between 0 and 1, of the equation 
n 
on ee 
n 


Accordingly, by a well known existence theorem for implicit functions,® there 
is a neighborhood 4 of « such that for each point z of A the equation 


o(0,2z)=0 
has a solution @ interior to the circle detined by the inequality 
@—4 <4. 


But for z +, (0, z) vanishes in the same points as @&(@, z), while /'(@, a) 
is identically zero. Our theorem is thus established. 

This method of demonstration also shows that, if we except the case 
where f(z) has the form az + 8, each branch of the implicit function @ is 
analytic when z is in the neighborhood of 7, and approaches one of the 


1 
(n— 1)th roots of = when z approaches a. 


2. Restrictions on Theorem I, and a Related Theorem. It 
should be noted carefully that we have proved only the existence of a region A 
in which the above theorem holds,—nothing has been said as to its size. In 
fact the determination of such a region for a given function and a given point 
a is usually a difficult problem. We can easily construct elementary functions 
such that, vo matter how small b—a’ may he, when b has once heen fixed 
there is no value of 0 for which F(0, 4) vanishes, and again functions for which 
no value of 0 satisfying the equation F'(0, hb) = 0 lies inside the circle defined 
by the inequality 


The function e*, where 





*For a clear statement of this theorem and references see Encyclopddie der mathe- 
matischen Wissenschaften, II B1, p. 103. 
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is an example of the former sort, while an illustration under the latter case is 
furnished by the polynomial 

z—a - 

a ”) . 


where m is an integer greater than 3, and ¥ is the solution 


7= a(1 + ‘cot =) 
m 


of the equation 
(y ae | = f*. 


For this function, @ = y is the only solution of the equation F'(6, b) = 0, and 
we have 


A closely related question is involved in the discussion of the equation 
(5) F(9, 2%, 22) = f(%) —f(22) — (4 — 2) S' [22 + 8 (4 — %)] =9, 
where both z, and z, are variable. If f(a) + 0, wecan make use of an auxil- 
iary function defined by the equations 
F(6, 215 22) 


a ° , 
~ > - 
(21 — 2 


(0, 2%, 2;) = (4 — 8) f"(4). 


This function vanishes at the point (4, a, 2) and is analytic in the neighbor- 
hood of that point, while 


(6) (6, 2 2) = 


~_ s+» 
“1 ~ “2s 


é 
oe 0. 
ap P (4, a,a) + 


We can therefore state the following theorem : 

TuHeoreM II. Let a be an interior point of a region B throughout which 
Sf (2) is analytic, while f"(a) 40. Then there exists a neighborhood A of a 
such that for all values of z, and z, in A the equation 
(5) S(a) —S (a) — (4 — 2) Pla + 8 (a — &))] = 0 
has a solution @ satisfying the inequality 
@—43\ <4. 
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ia ie 
an ae et If f(a) = 0, the auxiliary function ® detined by (6) cannot be used, and 
aan | a in fact the ahove theorem may he false. For example, if f(z) = 3s", Ww here m is 
me oS an integer greater than 3, it is possible to find points 2,, z, in every neigh- 
a # 44 borhood of the origin such that all solutions of (5) satisfy the relation 
ot 2 oe ern 
L 1S } ae er 
at 44 @—4\/ 24. 


We have this result when 2, =e ™ z,, since the only solution of (5) is then 


we 
: ie ~om wae = 
: 
a 
2 
a 


— 
or 


a, l 


seedentt egies 
‘ . 


‘Bi 

Ea ae ie l—em 

Kia 

eG Es 3. Formulas Involving Integrals. As a preliminary to the de- 
ite 2 ° 
he velopment of theorems more general than the preceding, we shall now consider 
He 


* 


siete 
Mey 
Bowne eh) yerteons 


the equation 





eae a 


“ 


‘ . 7 "2 "a+t@(s—a 
a me a, (7) F, (0, z= | T( uw, 2) d(r, z) dw — f(a, z) / d(i, z) dw = 0, 


which is analogous to the formula of Bonnet’s theorem.* The functions 
T(v, 2), (1, z) will be taken as analytic functions of both w and z for all 
values of those variables satisfying inequalities 


be Na 


were 









ety eure 


sm ai <A. wo — @ < k, 


while f(a, z) will be supposed not identically zero in the region defined by 
the former of the above inequalities. 
For the functions f(1, z), $(7, z) we have developments 


f(r, z) = > Sa.(07 —a)*(z—a)’, 


u=z0 v=0 


d(w, z) -> >> B,,(w —a)*(z— a)’. 
v=0 


ast) 


= » ee - me 
Sra cet eit waitin ap Ph Py tat ta 
. . 


By substituting these developments in (7), we obtain for F,(8, z) the ex- 
pansion 


(8) F (6, z) => P, (0) (2 — ay*), 


A=O0 











* Journal de Mathématiques, vol. 14 (1849), p. 249. 
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where 
A 1 A—m “ 
> —_ 
(9) Pyx41(9) > pti pm > a,,B8, ee + 1000B,,.—n—e | 


For any value @ = @,, and for any positive real number A’, there exists a 
positive real number // such that the series (8) converges throughout the 
region defined by the inequalities ‘ 

\6-—6,| < A, z—a\<H. 


The following theorem can now be established : 


TueoremM III. Let f(w, z) and $(w, z) be analytic functions of buth w 
and z throughout a region defined by inequalities 


jz—a| <h, lwe—a\<hk, 


and let r, be the least value of r Jor which the polynomial P,. (0) does not 
vanish identically. Then if T is any region in the 0-plane of which a root @, 
of the equation P, .,(@) = 0 is an interior point, there exists a neighborhood 
A of the point z = a such that the equation (7) is satisfied by a value of @ in 
T for every value of z in A. 

The auxiliary function to be used here is obvious. If 6, is not a multiple 
root of P, .,(@) = 0, so that 


0 

60 
we use the same existence theorem for implicit functions as in the proof of 
theorem I. In the case of a multiple root we avail ourselves of the more 
general existence theorem which may be found, for example, in Picard’s 
Traité d’analyse, vol. 2 (2d. edit.), p. 261. 

From theorem IIT we can deduce formulas which correspond to familiar 
theorems of mean value. In the one which follows, an analogue of Bonnet’s 
theorem, it is interesting to note that as z approaches «, @ approaches a point 
on the unit circle instead of a point in its interior as in theorems I, II, and V. 

Tueorem IV. Let f(z) and $(2) he analytic throughout a region B of 
which a is an interior point. Then there exists a neighborhood A of a such 
that for any value of z in A the equation 


(10) [ ST (w)o(w)dw - sa [* d(w)dw = 0 


P.41(9) 4 9, 
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has u solution @ satisfying the inequality 
@—1/ <1. 


Equation (10) is a special case of (7) where f and ¢ are functions of 
but one variable, so that 


a, =8,,=0, vw #9. 


By hypothesis 
J (a) = Go + 9, 


hence, if 8,9 is the first non-vanishing coetticient 8,5, we have 


Py 41(8) = SPEC — +1), 

Since @ = 1 is always a root of P, ., (@) = 0, theorem IV follows from theo- 
rem III, the special region 7 assigned in IV having been chosen in order to 
present a formula as closely analagous as possible to the ordinary Bonnet 
formula. 

The two formulas of the following theorem are usually referred to as the 
first and second theorems of mean value for integrals : 

THEOREM V. Let f(z) ond $(2) be analytic throughout a region B of 
which u isan interior point. Then for each of the equations 


(11) | Feyo (eye —s(« + 2 — “)) | ¢(0) aw -_ 


“a+O(.—a 


(12) [feneerae rea) | }(w)dw—f(z) $(w)dw=0, 


B a+6(z—a) 


there exists « neighborhood A of « such that for any value of z in A the cor- 
responding equation has a solution @ satisfying the inequality 


<4. 


To deduce this result for formula (11) we make the following substitu- 
tions in equation (7) : 


0-4 


S(w, 2) =| : $,(w)dw, 


$(w,z) = f{(wv). 
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Integration by parts now gives to (7) the form 
F,(4, 2) = [fice [dea] + [Ace b.(0y dw 
- [Ala + 0 -«)] Ao) | J duende 
= [Aco si(e) dw - fila + (2-2) [by(wydw = 0. 


We have only to drop subscripts to obtain equation (11), and the part of our 
theorem which relates to this equation is demonstrated if we can show that a 
root of P,,4: (8) = 0 lies in the region defined by the inequality 


|9-3| <3. 
The functions f(w) and ¢(w) will have developments 


I (w) = a,(w — a)’, 
(13) a 
o(w) = > 2 B,(w—a)’. 
v=0 
If y(w) is a constant, or if ¢(w) = 0, equation (11) is obviously an identity. 
In any other case let «,, and 8, be the first non-vanishing numbers of the re- 
spective sequences 





a, ag, ee ty 
Bos Bi, :- « 
We then have 
y= m +h, 
1 em 
Py 41(9) = anhr(;. +k+1 k+ i): 


and this last equation always has a real root 0, satisfying the inequality 
|0— 3) <4. 
Equation (12) may be derived from (7) by the substitution 
f(w, 2) =f(w) —f(2); 
$(w, z) = $(w). 
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As in the discussion of equation (11), we note that if /(w) is a constant, or if 
$(w) = 0, (12) is an identity. In any other case a,, and 8, have the same 
meaning as before, when f(w) and ¢(w) have the developments (13). We 
find 


A, =m+h, 





; _ ame a — oF+) 
Pyi(9) = eallrees , i 


so that in this case also 7, ,,(@) = Ohas a root between 0 and 1 on the axis 
of reals. 


Cauchy’s formula * 


S(z) —f(@) _S'{a + O(z — a) } 


Find ee SA Bint sd: |0—3| < 
$(z) —G(a) P'la+ O(z— «)]) sliaios 
and the various forms of Taylor’s series with remainder are special cases of 


the first formula of theorem V. In particular, theorem I is « corollary of 
this part of V. 


NORTHWESTERN UNIVERSITY, 
Evanston, ILL. 





* Calcul differential, p. 87. 


























NOTE ON REGULAR POLYGONS * 
By CHARLOTTE ANGas Scort 


It is a well-known fact that the circle described with its centre at any point 
P of a rectangular hyperbola and passing through the other extremity of the 
diameter through /? cuts the hyperbola again at the vertices of an equilateral 
triangle. The object of the present note is to show that the circles which cut 
a rectangular hyperbola at four of the vertices of a regular pentagon, hepta- 
gon, or nonagon can be described with equal facility. 

It is necessary, in the first place, to find a rectangular hyperbola that 
shall cut a given circle, x? + y* = 1, at four vertices of one of these polygons. 


1. The Heptagon. If a side of the heptagon subtends at the centre 
an angle 6, then 7@ = 27, hence sin 36 = — sin 40; that is, 
3 sin 6 — 4 sin® 6 = — 4 sin 6 cos 6 (2 cos? 6 — 1), 
from which, 
sin 8 }8 cos* @ + 4 cos? @— 4 cos @—1/ =0. 


If one vertex is at (1, 0) (see figure 1) the roots of this equation give the 
values of @ for the vertices. If x is written for cos 6, the solution sin 6 = 0, 
which refers to the vertex (1, 0) and gives therefore @ = 0, becomes 
x —1=0; and the cubic factor gives 

823 + 4x? — 4dr —1 = 0. 

Hence the roots of 


(x — 1) (82° + 42? —- 4x —1) = 0, 


that is, of 8at — 4273 — 82? + 382 +1=0, 
are the abscissw of the vertices numbered 7, 1 and 6, 2 and 5, 3 and 4, in 
figure 1. 


It is desired to find a rectangular hyperbola that shall meet the circle 
at the point 7, and at one point of every pair (1, 6), (2, 5), (3, 4). There 
are clearly 2°, that is, 8 such rectangular hyperbolas; any one will serve the 
purpose. 








* Presented to the American Mathematical Society, October 27, 1906. 
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The general equation of a rectangular hyperbola is 
ax? + 2hay — ay? + 2gx + 2fy+ec=0. 
To find where this meets the circle 
et+y?—1=0, 
eliminate 7. The two equations give, by elimination of the term ,/*, 


2ax? + 2hey + 2gx + 2fy+ec—a=O0, 
that is, 

Jaz? + 2Qgxa+c—-a=—y (2hx + 2); 
hence, since 


—-y—2z?—1, 
(x? — 1)(2hx + 2f)? + (2ax* + 2gx + e¢—a)?=0. 
The equation for the abscisse of the common points is therefore 
(4h? + 407) xt + (RAP + Bag) 2 + (4f? — 4h? + 4g? + 4a(e — a) (x? 
+) — 8Af + 4g(c —a){ x — 4f? + (ec —a)?=0, 
and this is to be the same as 
8a! — 4a? — 82° 4+ 37 4+1=0. 
Hence, equating coetlicients, 
4h? + 4a? 8A, (1) 
8hif+ ag =-- 4A, (2) 


4f? — 4h? 4 49? + 4a(c—a) =— BA, (3) 






—hf+4g(e—a) = Ba, (4) 
—4f? +(c-—a)y?= 2X. (5) 






The sum of these gives 









49> + 4g(c +07) + (e+ a)?=0, 





that is, 2q +c+a=(), 
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The sum of (2) and (4) gives 
4y(c+a)=—2X. 
Hence —%Wg=a+e, A= 2(a4c)?. 
With these values for g and A, (1), (5), (2) become 
h? = 4(a + c)? — ad, 
4f? = (c —a)* — 2(a +c)? =— (a + €)? — dae, 
2hf=a(a+ec) — 2(a +c)’. 
Hence, in virtue of the identity 2? x 4/? = (2hf)*, we have 
}4(a + ©)? — |} — (a + €)? — dacl = (a +c)? Ja — 2(a + ©) }?, 
Since 2(a + c) — a is a factor of each side of this equation, one solution 
is 2(a + c) — a= 0, that is, 
a= — 2c. 
This leads to 29 =—(a+c) =¢, 
A = 2(a +c)? = 2e, 
W@=4(a+c)?-@=4?-42=0, 
4f? = — (a+ c)?—4dac=—C + 82% = 7c’. 


Since c is a factor in all the coefficients, it is not zero; hence we may give 
it any convenient value, for example, —1. Then a= 2,h=0, 2g =—1, 
2f = + V7; choose for 2f the value + y7- 

The four possible values for a: c, and the two for f, account for the eight 
solutions of the problem. The other three values for a: c, given by the cubic 


}2(a +c) +al\)—(a+ c)? — dac! = (a+ c)*}2(a+c)—al, 


are irrational. Since only one solution out of the eight is needed, it is unnec- 
essary to attend to these. 
The rectangular hyperbola is 


247 —- 2yy —-x+Viy—-—1=0. 
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The center is at (+, + V7); when the origin is transferred to this point the 
equation becomes 


hence the semi-axis major is} ¥2. The curve lies as shown in figure 1; it 
passes through the vertices 1, 2, 4, 7. 


Since the foci of x? — y? = 3 are (+ $, 0), with the original axes they are 
(3, V7). (— 4,47); of these the first lies on the circle. The corre- 





Fic. 1. 


sponding directrices are x = + 4, that is, referred to the original axes, 


x=t++4=4,andx=—}4+4=9: the second of these passes through the 
centre of the circle. The distance between the foci is equal to the radius of 
the circle. 


2. The Nonagon. In the case of the nonagon, the angle @ ix deter- 
mined by 36 = % 7: hence cos 30 = — $, that is, 
4 cos? @— 3 cos 864+ 4=09, 
82> — 6x +1=0. 
The roots of this equation are the abscisse of the vertices 1, 4, 7 (see 
figure 2), hence of the pairs of vertices (1, 8), (4, 5), (7, 2). To bring in 
3 and 6, for which x = — 4, we multiply by 22 + 1; the equation becomes 


l6a4 + &e? — 127? -~dr+ 1 =0. 


By the same argument as in the case of the heptagon, the coefficients in 
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the equation of the desired rectangular hyperbola are determined by the five 
equations j 


4h? + 4a? = I6,, (1) 

8hf+ Bay = 8&2, (2) 

4f? — 40? + 49? + 4a(e — a) = — 122, (3) 
— 8hf+ 4g(c —a) =— 4), (4) 
—4f? + (c—a)? = r. (5) 


These, combined exactly as before, give 
49? + 49(c +a) + (c+ a)? =92, 
4g(c +a) = 4), 


hence 4g? — 5g(c + a) + (c+ a)? =0, 
that is, }4g— (+ 4){ }g— (e+ a) =0 
from which 4g=c+a or g=cr+a. 


Of these, 4g = c + a leads to the simplest results; we find 44 = (c + a)’, 
and then equations (1), (5), (2) become 


= (a+c)?—-@=c(2a+c), 

16f* = 4(c — a)? — (a +c)? = (a — 3c) (3a — Cc), 

4hf = (a+c)?—a(a+ c) =c(a+c). 
Hence, by means of h? x 16/7 = (4h/)?, 
we have 

c(2a + c)(a — 3c) (3a —c) = e(a + €)?. 

One solution of this equation is c = 0; the other three are irrational. We 
take therefore c= 0; then hk = 0, 4g = a, 16f? = 3a*, from which 
4f=+4 V3a; take 4f= + V3 a. 

(Since there are four pairs of vertices, the number of solutions to be ex- 
pected is 2‘, that is 16. All are accounted for; two values for g, four for 
a:c, two for f, give 2 x 4 ~ 2 = 16 solutions). 

Since a is a factor in all the coefticients, it is not zero; hence we may 
assign any convenient value, in this case 2. The desired rectangular hyper- 
bola is 
2a? yi a+ yd y=0. 
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The centre is at (— 4, 4 V3); when the origin is transferred to this point the 
equation becomes 
l 


alii ila! 





hence the semi-axis major is } ¥2. The curve lies as shown in figure 2; it 
passes through the vertices 2, 3, 5, 8, and through the centre of the circle. 
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Fig. 2. 


The radius of the circle to the vertex 3 is a diameter of the hyperbola; the 
distance between the foci is equal to the radius of the circle. 

The rectangular hyperbola thus obtained for the nonagon is of the same 
size as that obtained for the heptagon, but differently placed. 


ne ie — ates 





3. The Pentagon. In the case of the pentagon, 50 = 27, hence 
sin 26 = — sin 38, 
from which 


2 sin 0 cos 0 + 3 sin 0 — 4 sin? @=0; 
that is, 


+ a 
i 
ioe 
_ * 
aad: 
fp 
iis 
»* 
es Be 
2 


sin@=0 or 2cos 0+ 3 — 4 sin?@= 0. 
Hence 4 cos?@ + 2 cos8d—1=0, 
that is, 427+ 2zx7-—1=0. 


ne 





The roots of this quadratic are the abscisse of the pairs of vertices (1, 4), 
(2, 3) see figure 3. A rectangular hyperbola, with the axis of x as one axis 


of symmetry, can be passed through the four points. 


The equation of such 
a hyperbola is 


ax?— ay + 2ge+c=0; 
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the abscissve of the intersections with the circle are the roots of 
ax* + a(x?—1) + 2gxr+c=0, 
that is, of 2a? + 2gx — (a—c) =0. 
Hence 2a = 4, 2g = 2, a~—c = 1, and the hyperbola is therefore 
Qa? — dy? + We 4+1=0. 


The centre is at (— 4,9); when this point is taken as origin the equation 


becomes 
xt — ¥ + } = 0, 


hence the semi-axis major is $. The curve lies as shown in figure 3. The radius 
of the circle to the point (— 1, 0) is the minor axis of the hyperbola. 


14 


to 








Fic. 3. 


4. Conclusion. If we change the origin in all three cases to the centre 
of the hyperbola, and change the scale so that the hyperbola has the semi-axis 
major = 1; and moreover, in the case of the heptagon, interchange the axes 
of x and y; we arrive at the conclusion that the hyperbola x? — y? + 1 = Ois 
cut at four vertices of a regular pentagon by the circle (x — 1)? + y= 4; at 
four vertices of a regular heptagon by the circle 


iT s w 4 ~) = 8: 
(= - 5) + (¥+ a) = 8 


at four vertices of a regular nonagon by the circle 


(*~ ¥)'+ (v+ V5) = 















ed 


ee eee 








SS. OE OS TT Pe: 5a emmerete 


itn maa al 


? 









ER OO RR Pr 


Cr 


- —apenaeeionnyn es menylpnepeamnttinet 2 samaphiagorgs ta 


wg 








134 SCOTT 


_The relation of the different circles to the hyperbola is, however, more neatly 
expressed by the geometrical statements of the results : 


(1) A rectangular hyperbola is cut at four vertices of a regular penta- 
gon by a circle on the minor axis as a radius (figure 4). 
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(2) A rectangular hyperbola is cut at four vertices of a regular hepta- 
gon bya circle of radius equal to the distance between the foci, with its centre 
on one directrix and passing through the other focus (figure 5). 

(3) A rectangular hyperbola is cut at four vertices of a regular nona- 
gon by a circle which has for radius a diameter of the hyperbola of length 
equal to the distance between the foci (figure 6). 
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THE REVOLUTION OF A DARK PARTICLE ABOUT A LUMINOUS 
CENTER 


By Epwin BipweLL WILson 


1. Introduction. The simple problem of two bodies moving under 
their mutual gravitational attraction has doubtless long been considered as 
completely solved. Yet ever since Maxwell’s formulation of his theory of 
electromagnetic action as a motion or stress of some kind in the ether, there 
has been a remote possibility that the time might come for revising the dynam- 
ical equations of motion so as to take into account the action of the ether on 
the moving masses. In fact, recent investigators of the electronic theories of 
matter have paid a great deal of attention to the interaction of matter and 
ether. To adapt their work to the problem of the revolution of bodies about 
a fixed luminous center, which is approximately the case that obtains in the 
solar system, is by no means easy. For in the first place there are different 
theories of electrons and in the second place these theories lead one to try to 
replace the older and more familiar dynamical equations by newer and as 
yet not fully settled electrodynamical equations. Thus Abraham™* discusses 
in considerable detail the actions of the ether on a moving body; but his 
equations are those of the electromagnetic field and depend to some extent 
on his particular theory. 

Recently, however, Poynting, with that directness of insight and sharp 
mechanical instinct which characterize the English investigators of physical 
science, has taken up the problem of the effect of radiation on moving masses 
and has treated it from the point of view of the transmission of momentum by 
the ether. His results, although perhaps essentially electromagnetic in that they 
are derived from the consideration of waves in the ether, are nevertheless 
well adapted to direct mechanical formulation.t The effects which he describes 


* Theorie der Elektrizitdt, vol. 2. 

+ We do not wish to be interpreted as implying either that the pressure of radiation is 
necessarily connected with any electromagnetic theory or that Poynting was the first to con- 
sider the matter from other points of view. On the contrary the question is rather one of 
mechanics and thermo dynamics and had previously been treated as such by several authors. 
Historical references may be found in an article on ‘‘ The pressure due to radiation ~ by Nichols 
and Hull, Astrophysical Journal, vol. 17, June, 103 or The Smithsonian Report for 1903, pp. 
115-1388. Poynting merely reverts to this simpler treatment and applies it anew. 
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are three:* the direct pressure of light, the Doppler reception effect, and 
the Doppler emission effect. We will state what is meant by each. 

If a body which is under the rays of the sun be considered as darkf or 
essentially black, to use the technical term, the radiation from the sun is 
stopped and absorbed by the body. The result is an outward pressure due to 
the abstraction of momentum from the ether. As the electromagnetic radiation 
varies in its intensity according to the same law as the gravitational attraction, 
it produces a pressure which varies inversely as the square of the distance and 
hence introduces no change in the dynamical equations except an apparently 
diminished value for the gravitational action. It might seem at first sight 
that in estimating the pressure due to radiation the distance should be meas- 
ured from the radiating surface: but a moment’s consideration shows that a 
spherical body like the sun could be replaced by a concentric body radiating 
with a properly altered intensity, and hence that the distance should be meas- 
ured from the center. 

The two Doppler effects are much smaller than the direct pressure of light. 
To consider the reception effect suppose that the dark body were moving 
toward the sun or any other luminous center with a velocity v. Then a 
greater amount of momentum would be absorbed per unit time than if the body 
were at rest. In particular the amount would be greater in theratio v/ V, where 
V denotes the velocity of light. In case the body falls to the sun from an in- 
finite distance this ratio is about 2/1000, which may be considered as its maxi- 
mum value in the solar system: generally it is much smaller. If the body 
retreats, the pressure is lessened in the same ratio. Finally there is the emis- 
sion effect which arises from the fact that the dark body radiates the energy it 
receives. This foree, which is “frictional” in that it always opposes the 

*See his address to the Royal Society, Vature, vol. 75, pp. 90-93, November 22, 1906. 
Earlier than this, Poynting had printed an extensive article on ‘‘ Radiation in the solar system : 
its effect on temperature and its pressure on small bodies,” Transactions of the Royal Society, 
ser. A, vol. 202 (1903), pp. 525-552. In this article there occur some minor errors which 
the author corrects ina note at the end. Moreover the two Doppler effects are by no means 


clearly distinguished one from another. Finally in this earlier treatment, loc. cit., p. 549, 
the author comes upon a differential equation that he states is probably not integrable 
This equation corresponds to our equation (5). The integration which we give in terms 
of Bessel functions not only furnishes a desirable complete determination of the problem, but 
enables us to see certain other general characteristics of the motion than are obtainable from 
a less exact integration. 

+ This is approximately true for the small particles in the solar;system. As far as we 
know them, their albedo is very low. 
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motion, is stated by Poynting to be one-third of v/V of the direct pressure of 
light.* This does not make the emission effect merely one-third as great as 
the reception effect: for the velocity which enters into the determination of 
the reception effect is the component velocity of the motion resolved along the 
line joining the particle and the sun, whereas in the emission effect it is the 
total velocity in path that enters into the formula. 


2. The differential equations of motion and their integrals. 
These statements of Poynting’s may be incorporated in the dynamical equations 
in a simple manner if polar coordinates with the sun as center be assumed for 
reference. Let — Cr? and A‘/r* respectively be the accelerations due to the 
gravitational attraction and the direct pressure of radiation. There will be the 
accelerations 

KW 1 dr 1 ff 1 dr 
-s FP and a oe 
r? Vv dt 3? V dt 
along the radius vector due respectively to the reception and emission effects. 
Further there will be the single acceleration 


lf 1 rdé 
3r V dt 


perpendicular to the radius. The equations of motion then become tf 


a dée —~K 4 Hh dr 

(1) -(G)=- 3 Vr dt’ 
dé t rdo 
(2) r : 5 a (" 1) 34 Pa “dt 


The simplicity of these equations is sitedied due to the fortunate linear 


*As a matter of fact Poynting states, as is well known, that higher powers of v/V 
enter into the expression; but neglecting them does not introduce, on the assumptions 
made above, an error of more than one-fifth of one per cent in a quantity already small; 
and generally the error is much less than this. To introduce the exact expression for the 
emission effect (which is evaluated on assumptions that are not perfectly well established) 
would unnecessarily complicate the equations of motion. 

+ Provided we disregard, as we may, any changes in the mass of the particle which 
might arise from its motion through the electromagnetic field. On the electronic theories 
there is such an effect, but it remains very small until the velocity of the moving body is- 


within a few percent of the velocity of light. 
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dependence of the Doppler effects on the velocity of the particle. Equation (2) 
is immediately integrable into 


, a0 iM 
(3) Mat 


dt 3 0. 





To find the orbit of the particle we may eliminate the time between this equa- 
tion and (1). Meanwhile in order to reduce the resulting equation to a 
simpler form, it will be well to introduce the change of variables 


1 3VA 
= — —+- —9=¢. 
(4) u=-, $ 
The differential equation of the orbit then takes the particularly simple form 
5 @u 3 du oa 9? C-HK L 
(5) mm -2H**"*3-F "sR 


Two things may be noted at this point. First, if A’ be allowed to ap- 
proach zero, ¢ becomes infinite in such a manner that Ad remains finite and 
dé = —d0@. Hence (5) approaches, as it should, the form 


au 


102 + uz const., 


ordinarily obtained in the theory of gravitation for the revolution of a planet. 
Secondly it is noteworthy that (5) is, except for the presence of the second 
member, an equation of the Bessel type. 

To bring the differential equation of the orbit into shape for easier com- 
parison with the treatises on Bessel functions the change of variable v = z¢* 
will be made. Then (5) takes the final form: 


; d?z 1 dz 2? L 
(6) et ae +(1-3)2=3- 
dg? odd ¢* ¢' 
If the second member were zero, the general solution would be 


z= AS, (¢) T B Y,(¢) 


in terms of Bessel functions J and Neumann functions ).* It remains to find 
a particular solution of the complete equation to add to this. These particular 


*The notations throughout this article will conform to those of the Treatise on Bessel 
Functions by Gray and Mathews (1895). 
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solutions, rarely given in the books, may readily be obtained by the method 
of variation of constants. In this case such a solution is , 


Hence the equation of the orbit may be written as 








u=g | Ah($) + BY) + L[ V9) | oa ag — Jy) | 4) ag] |, 
where the constants A and 2 may be chosen so that the lower limit of inte- 
gration may take any desired value. 

As the Doppler etfects are very small, the motion of the particle will be 
nearly elliptic, parabolic, or hyperbolic. To bring out more clearly the de- 
partures from these forms we may expand the solution (7) by means of the 
asymptotic expressions for the Bessel and Neumann functions. These series 
will semi-converge with extreme rapidity because of the large value of the 
argument ¢ (we shall have no occasion to consider cases in which ¢ is mat- 
erially less than 10° as may be seen by reference to p. 148). For such cases as 
these the function Y is inconvenient: it contains an additive multiple of J. 
Let us write * 

VCP) = (log 2 — 7) In(b) + Zn (¢)- 
Then the function 7 satisfies the fundamental relations which J and Y satisfy 
and which are the reduction formulas} especially useful in the manipulation L 4 
of these functions. We have then the asymptotic formulas 


m\!* “ (2u + 1)7 } ‘ ((2n + 1)r | 
(Y) J,(¢) = 1 cos | — or + a) sin | d 


| 
4 4 j 





26\1 2, ((2n+1)7 .) _ ((2n+1)r_ 4) 
—_— 4 4 ‘Os <. ——c—“—,—-, — | —- Ps < -—- -—-—-— — a) : 
(=) Zn() yon 4 e In { 4 % | 
where 
4n? — dn? — { 1 

P =i — (= 5 reaps és; + higher even powers of $° 

4n? — 1 (4n? — 1) (40? — 9) (4a? Ha 25) higher odd powers of : 

C= ag ~ 81 8G)F aos wen ¢ 


*See Gray and Mathews, p. 40, formula 95. 
+ These are: Z', = (n/¢)Zn — Zn41 =Zu—1— (n/O) Zu, Zn41— (2n/O)Z, + Z,-1=9, 
Jne1 Zn —InZn41=0, and (¢"Z,)' = 9"Z,—-1. See Gray and Mathews, p. 14 and p. 16. 
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Or if we set (2n + 1) 7 4—$= , We may write 


; -_ mn 2 
(7) "J,() a ma + n " } sin vi, -_ (n oar - #) cos Vn end ese 
26\!2, 


with the relations cos y,_,;=sin y,, sin ¥,_; =— cos fy. The terms 
omitted from (8) are of the order 10- when ¢ is of the order 10°. 
In terms of the functions J and Z the equation of the orbit becomes 











Il 


1 t cos ¥, + (? sd ld siny, +--+, 





— sin yy, > Sa 


v= G1 CAL) + DA) + L[ A) [Da — 1) [AYPag] |, 


where C’and D are constants of integration slightly different from A and B. 
This form of the equation is not yet in the shape most convenient for computation. 
For the integrals in the bracket change their value by a hundred percent when 
¢ runs through an interval of 27. It may be noted, however, that the maxi- 
mum value of the integrals, when extended over a small interval in ¢, is very 
small owing to the fact that the denominator remains of the order 10" and the 
numerator of order 10%. An integration by parts will enable us to ac- 
centuate this. We have 

HEB) ip = [PHD ag PHD) 5 5 (PHD) ag 

d ¢° J ¢ 

JS3() ¢* J3($) 
Ae) 6/40) ap 


I) J3() S,($) 7 J5(p) 
[AG = GP + GP + 8 


with a similar result for the integral of 7,(¢)/¢*. Combining, we have for 


the bracket 
4s 
+ higher even powers of 


2(b) 49) . 
Z 4 — J, _ < ; : 
Wo) [3 ” (6) vats > ¢° ¢ 


There appears to be no necessity of keeping terms of the order 10-*. Hence 





Hence 














fee, 





(10) u = g#! | CAG) + Déal$) + L( 5 - #5) } | 


cu A RM 
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may be taken as the approximate equation of the orbit for future compu- 
tations. For the derivative we have, from (10), the expression 


(11) w= $1 C1(b) + DA() — L at 


where the small quantity 48 Z/¢’ has been nagieoted within the brace. 
If the initial values uo, u'o, and do of u, u’, and ¢ be substituted in (10) 
and (11), the constants C and D may be obtained. 


eo b{le-H- Daw -[ 
= 5 (RA(b) - SZ,($0)), 
D= al - [ — u(s - 5) vue) + [ + =] ($0) | 


= 5 (- Bs) + Sh(0))- 


In equations (8) let y, = ~, = y. Then we have, if (2n + 1)7/4— $ = Yo, 





Zs (40) } 


(“F)* Si($o)=sin Yo — 54- C08 Yo — hae sinyot---, 1 iA 


Ce)" Z1($o) = COS Yu + s sin Yo — er sa Yotrcr:. 


Hence 


Zu($0) (6) = ) is Yo cos Y + COS Yo sin v5 + sin Yo cosy 3 


and 


————~ 


Si ($0) Z2(%) = (a 3) {- sin Wo sin Y + sin Yo C08 Hot + COS Po siny - ora 


- cos y, cos p * ae ~ sin Yo sin ¥ aah ~ sin Yasin ¥or9g. 2 | 
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4 q Hence 


) 1\12 
RZ,($))42(¢) — RA (bo) 22(¢) = F (33) [ cos (yo —v) 


. 15 3 105 
— sin(vs — ¥) (55 - 3g, + cos( Wo — »(a3- 123g? * isa) |: 


In a similar manner 


— SZ -(b0)42(¢) + Sdo( bo) Z2(¢) 
12 
= 8(55) [singe v) + £08 CY — 9) (5 ~ 4.) 


pie ; 5 

+ , 225 105 105 
e4 +h 77 (sae ~ Tax? aa) | 
BF. 

| ‘s Now let —w=&. Then ¢ — ¢, = Eand = ¢g (1 + E/bo). Hence, by 
: 


expanding into series for two terms and collecting, we have 
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fase =| + 18(5- 3): 


It should next be noted that the coefticients R and S consist of the sum of 
terms of different orders of magnitude. In fact u, — 1/6, may be considered 
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Hier bo as of the second order whereas 12/¢! is of the fourth. So also in S we have 
1 bp a u, of the second order (unless the orbit differs only infinitesimally from a 
its circle, wv’ and wu are of the same order) whereas 2 ¢° is of the third. If we 


rearrange the expression for u on this basis and at the same time expand 
($/$o)*? we get as the final form. 


c= [(u—Z)oore + wine + 2] 
0 


(12) “* aE Gs - 31) sin & — a E+ 5 Mag sing - 3(% — aE cos € | 


+3[- = U4 ma sin’ $£ — “Pu Ecos & — (ua — ae sing 
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and hence, on differentiating, 


v= [ cos — — (wo _ =) sin é| 


(13) + zl u (sin E+ Ecos é) _ 9 (% _ 4.) ésiné — sin? 4] 


+ 1 F6L F 3 ss 69 L i 
“~~? : = 0 Ss — = en 1 
gL hi, 8 8 d, 


’ ; (« 4 =) E cos £ + 3(« cos & — (% ~ =) sin é) ap 


where the first, second, and third brackets denote respectively the first, second, 
and third approximations. The integration of the problem, both in exact 


and in approximate form, may therefore be considered as accomplished. 


2 F 
E+ my & sin & 


3. Discussion of the characteristics of the motion. In the first 
place let it be assumed that at the start the areal velocity is positive. Equation 
(3) shows that as the motion goes on and @ increases from its initial value 0, 
the areal velocity will steadily decrease and hence by equation (4) the argu- 
ment ¢ will decrease. If we allow ¢ to approach zero in the complete in- 
tegral (7), it appears that u is becoming infinite: for the terms 


>) dd 





$*J,(b), $7 )3(¢), 14(6) [ 


approach finite limits and the integral* 


J,(¢) 
d 
[Ape 


diverges for ¢ = 0. Hence the particle will reach the attracting center when 
@=0. It should be noted, however, that in an actual case the motion will be 
stopped before the particle reaches the center by the impact of the particle on 
the surface of the attracting mass. For another reason it is also valueless 
to investigate what happens to the orbit for very small values of ¢, that is, 
when the particle is near the attracting center. We have assumed in setting 
up our equations that the velocity of the particle did not become very great. 








* This follows immediately on substituting the principal values of the functions ob- 
tainable from the developments in series. See Gray and Mathews, p. 11 and p. 14. 
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Hence the particle should not be allowed to approach the center too closely. 
Hence 3 VA/X is an upper limit for @, and we may say that the number of rev- 
olutions that the particle can make is less than 

14 V = pate 

ang v= oR 


This number for any given particle depends solely on the initial areal velocity 
and not on the velocity of projection or the initial distance from the center, 
except in so far as they conjointly determine the areal velocity. 

The approximate solution (12) will serve to discuss the motion for a few 
dozen revolutions from the initial position ¢ = $9. If it be assumed that the 
polar axis is so chosen that u, = 0, that is, so that the particle starts from an 
apse, then the first approximation may be put in the form 


L 
i (1+ © 008 €), oz ye. 
$i, 

This is a conic section as might have been expected. In particular it is that 
conic section which the particle would described under the action of the gravi- 
tational attraction diminished by the action of the direct pressure of light and 
free from the Doppler effects. This conic I shall call the osculating conic 
of the orbit. 

To find the amount by which the particle falls in from the apse toward 
the center during the first complete vevolution it is merely necessary to sub- 


tract the values of the radius vector as given by (12) for & = 0 and € = — 2m. 
Hence to the second approximation we have 


1 1 3 1 1 3 
Anu, = Uy N Ise (2 _ 5) e). Ar, = —- No Nive (2 _ 5} e). 


This quantity A ro is always negative, showing an approach to the center in all 
cases,” and varies between the limits 7)/4.V and 2r,/V as the osculating conic 
varies from a parabola to a circle. It is to be noticed that, if the eccentricity 
is small, the particle falls in toward the center by an amount about twice as 


ul 


“z=: 





* We have assumed tacitly that 1 > 0, that is, that the attraction is greater than the 
pressure of light. 
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great as that required to carry it to the center in NV revolutions if it moved in 
by the same amount during each revolution. This need not be perplexing, 
for according to the formula the particle is always falling in by a certain per- 
centage of the distance remaining (not of the original total distance). It there- 
fore appears necessary that the particle should fall in by an amount which at ; 
any rate for orbits of small eccentricity is decreasing from revolution to revo- - 
lution. In fact if the increment of the radius vector be found to another ap- 
proximation, it appears that such is the case.* 

As the eccentricity of the osculating ellipse seems to be of importance, it | 
appears worth while to calculate the change of eccentricity during a revolution. 
Formula (13) shows that if the particle starts at an apse, it will be at an apse 
after a complete revolution provided we disregard the effect of the last ap- 
proximation. This fact will materially simplify the discussion of the eccen- 
tricity. It enables us to write the equation of the osculating ellipse after one 
revolution as 


L L 
w= (totus — aes) om E+ aes 


On calculating the eccentricity of this ellipse by the usual formula we find 


lh =— ae t 
(15) Ae 9N ; 


Hence it appears that the ellipse is always losing eccentricity. To speak r 
roughly one may say that the orbit is losing its eccentricity from about two to , 
about fifteen times as fast as it loses its perihelion distance according as the 
orbit is approximately circular or approximately parabolic. 

Equation (13) shows that although there is not much change in the angu- 
lar position of the apse, there is some. Let us assume that the apse advances 
by the amount ¢, so that if & = 0 is a root of uw’, then & = — 27 —e is the 
next but one. It has been seen that ¢ is small of the order 1/.V?. The second i 
term of the first approximation gives a term ine. All the terms of the second 


*If I understand Poynting aright he states in Nature, loc. cit. p. 92, second column, { 
that the particle will fall in farther during the succeeding revolutious than during the pre- ‘’ 
ceding ones. As Poynting very seldom makes a mistake, I have given the matter a great 
deal of attention but fail to find any error in my work. Moreover Poynting, in his earlier 
paper, obtains a formula which is interpretable in much the same way as mine. Unfortu- 
nately I caunot find that he has printed any mathematical treatment of the subject since 
then. 
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approximation are too small to contribute anything. The first and fifth terms 
of the third approximation contribute quantities independent of «. The equa- 
tion for ¢ therefore is 


( L l E ie uC ,) 27] 
— —. = = » * > —. 82 . 
. 5) $3 L di 5 i 


l 1 9 3 360 1 9 3 tte 
— =(: > x¢) = $2 (: - x) | gre S, 


and there is a slight advance of the apse in the direction of motion. If the 
eccentricity is very small, that is, small of an order approaching 1/V? the apse 


Hence 





may advance over a considerable angle during a revolution. * 

So much for the discussion of what happens in the vicinity of the start. 
We are now in a position to take up more general considerations. From the 
behavior of the apse investigated above, it should appear that as ¢ decreases 
from its initial value ¢9, the apses draw further and further apart; and con- 
versely if we follow the motion backward allowing ¢ to increase, the apses will 
follow each other at intervals diminishing toward the limit 7. Now this is 
what we should obtain on general considerations applied to formula (11). 
For the apses correspond to roots of u’, that is, to roots of the expression 


CMG) + DA() — 2. 


Here Cand D are definite constants. To be sure they are very small in the 
cases we have been considering, in fact of the order ¢;**. Nevertheless for 
values of ¢ sufficiently large, the third term, — 2,/¢°, may be neglected rela- 
tively to the first two. Then the roots of u’ may be investigated by the usual 
methods employed in discussing the roots of Bessel functions,t and we shall 
thereby obtain precisely the result concerning the advance of the apse that we 
have already obtained from our computations. As the value of wu’ becomes in- 
finite when ¢ = 0, it is clear that wu’ must have a definite smallest root corre- 
sponding to the last apse which the orbit possesses, after which it coils in 
around the attracting center with an essentially spiral form. It is not impos- 





*It is, however, by no means certain that the orbit would keep its approximately. con- 
ical nature for such small values of e: it might already have become spiral. See infra. 
+ Gray and Mathews, chap. 5. 
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sible to get a very rough estimate of where this last root occurs. For there 
can be no roots for values of ¢ so small that the third term above is larger 
than the other two. If we compare merely the magnitudes of ¢, ¢), C’, and D, 
we find that the condition is fulfilled when ¢ = ¢3°. To get any more ac- 
curate estimation seems not to be worth while. 

Similar considerations may be applied to formula (10) to discuss the roots 
of u. If the test suggested above be used, it appears that the roots of u are 
likely to be found for values of ¢ of the same order of magnitude as the initial 
value ¢). That this is correct follows from the approximate formula for uw. 
Here, if the eccentricity ¢ is greater then unity, there are roots of « in the 
immediate neighborhood of ¢ . If e is considerably less than unity, the 
motion of the particle may be followed back through a very large number of 
revolutions before the eccentricity shall have increased to unity. Nevertheless 
the rate of increase of e as given by (15) is tolerably fast, and it would be 
safe* to say that the first root of « would appear before ¢ had increased to 
many times ¢,. Now whenever u has a root, the particle has retreated to in- 
finity. Hence it follows that not only is the number of revolutions which a 
particle can make inward toward the center a finite number J, but the total 
number of revolutions which a particle can make in infinite time about an at- 
tracting luminous center is finite. This result is so different from that of the 
ordinary theory of attraction as to be very striking. In fact, it could hardly 
be credited if the analysis had not assured its correctness. For values of ¢ 
greater than the first root of u the physical interpretation is simply that the 
particle is describing a hyperbolic type of path. As ¢ gets larger and larger 
the roots of u become separated by intervals equal more and more nearly to 
am. In other words the path is becoming more and more nearly straight. The 
analytical deductions, therefore, are seen to fulfil in every respect the con- 
ditions anticipated by the physical intuition. 





*In fact if we express C and D in terms of e, and if we assume that the particle 
starts at an apse where ¢, has such a value that Z.(¢,) =0, then J\(¢,) =0 and the equa- 
tion for roots of u may be written roughly in the form 


le? 
‘a... L \t? 
1 Lely) =, to) =(55) - 


Hence we may look for roots of u in the neighborhood of 


fcc. oo 
7 = oy, & =o? “dy 
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In closing we may take the time to obtain numerically some of the 
quantities which we have been discussing. The case of particular interest is 
that of a particle revolving aboutthe sun. The following data may be assumed : 


velocity of light, 1 = 3 x 10” 

distance of sun from earth = 493 x V 

acceleration of earth in orbit = 59 x 10° 

pressure of sunlight at earth’s distance = 59 x 10°. 


All these quantities are in the C. G. S. system. The constant C in equation 
(1) has the value 


C=59 x 10° x 4937? V2? = 1-29 x 10”, 


The constant A may be determined in terms of the radius a and the density p 
of the particle as 


Kh: C = 59 x 10% wa? : 4rpa®59 K 10°, A= 3pta'lotCc. 


If the particle has a density about equal to that of the earth and a radius of 
about 0-O0lem, then A = 1-76 x 10%. It should be noted that A’ is increased 
by a decrease in the radius or density, and that the value of A’ here obtained 
is certainly small for such particles as meteoric dust finely divided. Suppose 
that the initial double areal velocity is 45 x 10'* which is roughly that of the 
earth. Then if a particle be started at an apse at the earth’s distance from the 
sun, 
bd = 2-3 XK 10', N= 5-66 & 10%, e= 0-052. 


Hence the distance of fall from an apse will be 80-6km during the first 
revolution. The change of eccentricity will be about one millionth of its value. 
The apse will advance about one half of one millionth of a second of arc. It 
is very doubtful whether the particle could make so many as 10° revolutions 
altogether. If the particle had a radius of one tenth the amount above as- 
sumed, the fall would be ten times as great and so would the diminution of the 
eccentricity : the rotation of the apse, however, would be one hundred times 
as great as it was before. 


YALE UNIVERSITY, 
New Haven, Conn. 
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